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ABSTRACT 

We present //ST/NICMOS [Fell] 1.644 /xm, Pao (1.87 /xm) and continuum images of the starburst 
galaxies M82 and NGC 253 at an unprecedented spatial resolution. In both galaxies we detect [Fen] 
compact sources superimposed on a diffuse background in the disk of the galaxies together with a 
component above and below the plane of the galaxy. The radio and [Fe n] emissions perpendicular to 
the disk of M82 show a remarkable similarity to each other, suggesting that both emissions originate 
in shocks from supernova explosions. We find a spatial correspondence between bright compact [Fe n] 
emitting regions and the location of radio supernova remnants (SNR) for approximately 30 — 50% of 
radio SNRs in M82 and NGC 253. This lack of a one-to-one correspondence, more than being indicative 
of a different origin for the radio and [Fen] emission in starbursts, suggests two populations of SNRs: 
an older population (<^ 10 4 yr) traced by the [Fen] emission and a younger population (a few hundred 
years old) traced by the radio SNRs. We therefore conclude that the [Fen] emission in starburst galaxies 
provides a good estimate of the supernova activity. Using our newly determined [Fe n] luminosities 
(corrected for extinction) of M82 and NGC 253, we reevaluate the calibration of the supernova rate in 
terms of the [Fen] luminosity for starburst galaxies. 

Subject headings: Galaxies: individual: M82, NGC 253 - galaxies: nuclei - galaxies: photometry - 
infrared: galaxies 



1. INTRODUCTION 

Compact radio sources in nearby starburst galaxies of- 
ten appear to be supernova remnants (SNR). SNR in 
our Galaxy and the Large Magellanic Cloud (LMC) are 
sources of bright near-infrared [Fe II] line emission with 
luminosities ranging from 10 26 W to 3 x 10 29 W (for the 
[Fe n] 1.644 /im line, e.g., Oliva, Moorwood, & Danziger 
1989; Burton & Spyromilio 1993). This behavior suggests 
a connection between supernova activity and the [Fen] 
emission in starbursts. This is supported by near-infrared 
spectroscopy (Moorwood & Oliva 1988; Greenhouse et al. 
1991; Vanzi & Rieke 1997), by high-resolution [Fen] imag- 
ing studies (Forbes et al. 1993; Greenhouse et al. 1997; 
Alonso-Herrero et al. 2000), and by the tight correlation 
between the radio 6cm and the [Fen] emissions (Forbes 
& Ward 1993). Further determination of the relation of 
the infrared [Fe n] emission with the supernova activity in 
nearby starforming galaxies would be helpful both for the 
understanding of these galaxies and for using the near- 
infrared [Fe n] luminosity as estimator of the supernova 
rate in more distant galaxies where individual radio su- 
pernovae cannot be identified. 

There are indications of some complexities in this re- 
lation. For example, the bright [Fe n] emission may not 
be spatially coincident with the radio SNR. Lumsden & 
Puxlcy (1995) conducted a spectroscopic survey for [Fen] 
emission from radio SNRs in M33 and found that in some 
cases the [Fen] emission does not peak at the position of 
the radio SNR, as well as that some radio SNR show no 
[Fe n] emission. Recently, Morel, Doyon, & St-Louis (2002) 



have carried out an [Fe ii] imaging survey of optically se- 
lected SNR in M33 and found a similar result. Both works 
have concluded that the [Fe n] emission from SNR may be 
produced over a time scale of ~ 10 4 years. Greenhouse 
et al. (1997) found that none of the brightest sources of 
[Fen] emission in M82 were coincident with the position 
of the radio SNRs. 

M82 and NGC 253 are the nearest starburst galaxies 
(distances d = 3.2 Mpc and d — 2.5 Mpc respectively). 
Both galaxies show a large number of radio compact 
sources (e.g., ~ 50 in M82 Kronberg et al. 1985 and Huang 
et al. 1994; ~ 60 in NGC 253 Turner & Ho 1985 and Ulves- 
tad & Antonucci 1991, 1997 and references therein). Their 
proximity makes these two galaxies ideal to study the re- 
lation between supernova activity and the [Fen] emission 
in starburst galaxies. Greenhouse et al. (1997) presented 
a Fabry-Perot [Fe n] 1.644 /xm map for M82 with a spatial 
resolution of ~ 1.3". Their main conclusions were that the 
compact [Fe n] sources contribute only ~ 14% of the total 
[Fen] emission in M82 and that they trace a population 
of SNR older than the radio SNR. Their work is some- 
what inconclusive because of the limited angular resolution 
and also because strong extinction could hide [Fe n] emis- 
sion from some of the radio SNRs. This paper presents 
fiST/NICMOS images of M82 and NGC 253. The high 
spatial resolution of the NICMOS [Fen] 1.644 /im images 
(FWMH = 0.35"; this corresponds to 5pc and 4pc for 
M82 and NGC 253, respectively) allows a detailed study 
of the location of the bright sources of [Fen] and of the 
role of supernova remnants in this line emission. 
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The [Fe n] emission in M82 and NGC 253 



2. NICMOS OBSERVATIONS 
2.1. Data Reduction 

HST/NLCMOS observations were obtained on April 7, 
8, 12 and 22 1998 (for M82) and on August 13 and 20 1998 
(for NGC 253) using the NIC2 and NIC3 cameras. Table 1 
lists details of the observations. The observational strategy 
consisted of taking images in a spiral dither with a 5.5 pixel 
spacing with four positions for a given field. The plate 
scales for the NIC2 and NIC3 cameras are 0.076" pixel -1 
and 0.20" pixel -1 , respectively. The FWHM of the point 
sources in the NIC3 F166N continuum images is 0.35". 

The reduction of the NICMOS images used routines 
from the package NicRed (McLeod 1997). A master dark 
image was produced by combining between 10 and 20 
darks for a given sample sequence after the subtraction 
of the first readout. The darks used were part of other 
proposals observed close in time, and the flatfield images 
were on-orbit data. The data reduction was performed us- 
ing the following steps: subtraction of the first readout, 
dark current subtraction on a readout basis, correction for 
linearity and cosmic ray rejection (using fullfit), and flat- 
fielding. The individual dithered galaxy images were reg- 
istered to a common position using fractional pixel offsets 
and cubic spline interpolation, and combined to produce 
the final image of each field. 

2.2. Flux- calibration of the images 

Prior to the flux calibration of the images, the back- 
ground needs to be subtracted (only for wavelengths long- 
wards of 2/im). Due to the large projected size of both 
galaxies, the field of view of the NIC2 images is not large 
enough to enable measurements of the background on 
empty corners of the images, so we used background mea- 
surements taken during the Servicing Mission Observatory 
Verification (SMOV) program for the filter NIC2 F222M. 
The F215N images cover a larger field of view, and an 
estimate of the thermal background was obtained from 
the corners of the images. The flux calibration of the 
broad-band, on-line and off-line images was performed us- 
ing the conversion factors (from ADU/s to mJy) based on 
measurements of the standard star P330-E during SMOV 
(Marcia Rieke 1999 private communication). 

Ideally one would like to measure the continuum on both 
sides of the emission line to estimate the continuum at the 
emission line. However, NICMOS only provides contin- 
uum bands to the red of the emission line. In cases of 
objects such as M82 and NGC 253, where the extinction 
to the nuclear regions is very high (even at infrared wave- 
lengths), a straight subtraction of the longer wavelength 
off-line image may result in an over correction of the con- 
tinuum at the wavelength of the emission line. To esti- 
mate the continuum for the NIC3 F164N images, we fitted 
the continuum on a pixel- by-pixel basis between 1.66 /im 
and 2.15 /im by linear regression of the flux as a function 
of the wavelength for the NIC3 F166N and NIC3 F215N 
line-free images. The extrapolated continuum was then 
subtracted from the NIC3 F164N image to produce the fi- 
nal continuum-subtracted [Fe n] 1.644 /im line image. This 
procedure does not correct for extinction, but it rather en- 
sures the subtraction of the correct obscured continuum at 
each position of the galaxy. For the NIC2 F187N image 
(Paa), since we did not have other narrow-band images 



taken with camera 2, we used the NIC2 F190N continuum 
image. Fortunately, the continuum need not be subtracted 
as accurately for Paa as for the [Fe II] emission line because 
the Paa emission is much brighter. 

The fully-reduced images of each individual field were 
combined to produce the mosaics. The images presented in 
this paper were taken during different orbits. This resulted 
in slightly different orientations (from P. A. = —100° to 
-105° for M82 and from P.A. = 57° to 68° for NGC 253). 
Before combining the images into mosaics, they were ro- 
tated by linear interpolation to north up, east to the left. 
This provided the same orientation for all the images, and 
more importantly allows comparisons with previously pub- 
lished data. 

To construct extinction maps to correct the [Fe n] emis- 
sion, and [Fell] /Paa maps, the NIC2 F187N, F190N, 
F160W and F222M images were demagnified to match the 
plate scale of the NIC3 images and filtered with a Gaus- 
sian with a = 0.65 pixel = 0.13" to match the point spread 
function (PSF) of the NIC3 images. 

3. THE [FEII]1.644/tM AND PAa MORPHOLOGIES 
3.1. M82 

Figure la and lb are mosaics (after rotation) of 
the continuum at 1.66 and the continuum-subtracted 
[Fe ii] 1.644 /im line emission, respectively. Both images are 
displayed on a logarithmic scale to stress the contrast be- 
tween the diffuse extended continuum and the compact 
sources. The [Fen] 1.644 /xm line map has been smoothed 
with a Gaussian filter with a = 0.1". The line emission 
image shows some residuals which are due to both the 
under-sampling of the NIC3 camera PSF and the varia- 
tion in the shape of the PSF with wavelength. In addi- 
tion, the right hand side of the line image shows a slightly 
lower S/N because the F164N and F166N images at that 
position had different orientations. 

The [Fen] 1.644 /im line emission exhibits two distinct 
components: a number of relatively compact regions dis- 
tributed along the disk of the galaxy and an extended dif- 
fuse component in the form of filaments extending above 
and below the plane of the galaxy. The filaments seen in 
the [Fen] 1.644 /im emission line show a remarkable simi- 
larity with the chimneys at the base of the superwind de- 
tected in radio (Wills et al. 1999), which are interpreted by 
these authors as expelled material from the central region 
of M82. 

In Figures 2a and 2b we present contour maps of 
the 1.66 /im continuum emission and the [Fen] 1.644 /im 
line emission, respectively. Comparison of the NICMOS 
[Fen] 1.644 /im line map with that of Greenhouse et al. 
(1997) reveals that some of the sources seen in the lower 
resolution map of Greenhouse et al. (1997) break up into 
individual sources at higher spatial resolution. In par- 
ticular, their source Fel contains some 4 or 5 individual 
sources. Source Fe6 (the most south-west source in Green- 
house et al. 1997 map) appears as a point source at our 
resolution in both the NIC3 F164N and NIC3 F166N im- 
ages, and the signal to noise on it is compromised by the 
differing image orientations (see above). Therefore it is 
not surprising we do not detect line emission from it. The 
NICMOS images detect a number of compact sources in 
addition to those reported in Greenhouse et al. (1997). 
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Table 1 

Log of the HST/NICMOS observations. 
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Note. — Column (1): NICMOS camera. Column (2): Filter. Column (3): Corresponding ground- 
based broad- band filter, emission line or adjacent continuum. Column (4): Galaxy. Column (5): 
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In Figures 2a and 2 b we display the positions of 44 radio 
SNRs with peak flux densities above 1 mJy (list provided 
by Dr. Karen Wills, see more details in Section 5.1). The 
distribution of the radio SNRs and their relation with the 
[Fen] 1.644 /im emission will be discussed in Section 5. 

In Figures 3a and 3b we display the mosaics of the 
1.90 /im continuum (NIC2 F190N images) and the Paa 
line emission (continuum subtracted NIC2 F187N images) 
respectively. The Paa image is dominated by a ring of ion- 
ized gas with an approximate diameter of 400 pc. The re- 
semblance with the mid-infrared [Ne n] line emission map 
at 12.8 /im (Achtermann & Lacy 1995) is remarkable. The 
south side of this ring of ionized material appears to be in- 
complete (see Achtermann & Lacy 1995), and it also shows 
a "hole" of emission at the center of the galaxy, which ap- 
proximately coincides with the peak of the 2.2 /im surface 
brightness. 

It has been suggested that higher extinction in the south 
half of M82 may be causing the uneven distribution of ion- 
ized gas in the ring. However, the ring appears to be com- 
plete in [Fen] 1.644 /im line emission. If extinction were 
the cause, we should expect to see a similar morphology 
in the [Fe n] emission, which is more affected by extinction 
than the Paa emission. It is therefore tempting to suggest 
that the south part of the ring of ionized gas is actually 
broken, perhaps due to SN explosions. This hypothesis 
would also explain the excess of [Fe n] emission in this re- 
gion. The presence of the hole of gas emission at the center 
of M82 along with the anti-correlation between the mor- 
phology of the hydrogen recombination lines and the CO 
index has been invoked as evidence for an outwards prop- 
agating starburst (e.g., Satyapal et al. 1997). The diffuse 
component of Paa emission is also in the form of "fingers" 
or "filaments" above and below the plane of the galaxy. 
It is now known from ground-based Ha imaging that this 



emission along the minor axis of the galaxy extends over 
llkpc (Devine & Bally 1999). Figure lb and Figure 3b 
show how there is a very close spatial correspondence be- 
tween the morphologies of the extended [Fe n] 1 .644 /im and 
Paa emissions. 

3.2. NGC 253 

Figures 4a and 4b are the 1.66 /im continuum and the 
[Fe n] 1.644 /im line emission maps of NGC 253. As seen 
in M82, the [Fen] 1.644 /im emission is in the form of a 
number of compact regions along the disk of the galaxy 
superimposed on a diffuse emission component. Figures 5a 
and 5b are contour plots of these two emissions. We 
also show the positions of the 6 cm radio sources with 
Speak > lmJy from Ulvcstad & Antonucci (1991), ex- 
cluding those sources that are classified as H n regions in 
Ulvestad & Antonucci (1997). 

Ground-based [Fen] 1.644 /im and Br7 maps at lower 
spatial resolution (FWMH ~ 1.3") were previously ob- 
served by Forbes et al. (1993). They found that the central 
[Fe n] emission of NGC 253 is dominated by three peaks 
of emission. Our higher resolution NICMOS image shows 
a larger number of [Fe n] compact sources. For instance, 
source B (in the Forbes et al. 1993 notation) breaks up 
into three individual sources. In addition to the nucleus 
and sources A and B detected by Forbes et al. (1993), 
there are other sources along the disk of the galaxy, as 
well as diffuse emission perpendicular to the disk, more 
prominent in the south half of the galaxy. We believe that 
Forbes et al. (1993) did not detect the diffuse compo- 
nent above and below the disk of the galaxy because of 
the lower S/N of their [Fen] map. We also mark in these 
figures the position of the peak (the filled square) of the 
near-infrared emission (also coincident with the peak of 
the Paa emission) which was identified as the nucleus of 
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The [Fe n] emission in M82 and NGC 253 



Table 2 

Extinction to the nuclear regions of M82 using a foreground dust screen model. 
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Note. — Apertures are circular and centered at the peak of the 2.2 fim continuum brightness, except for sources El, Wl and W2 in Achtcr- 
mann & Lacy (1995) notation. The approximate locations of the Bra (or Pact) sources from the peak of the near-infrared emission are: El 
(+4.1", +0.5"), Wl (-4.6", -2.0") and W2 (-10.3", -4.5"). The errors quoted for the extinctions include uncertainties from the aperture 
size and centering, and flux calibration. 



Table 3 

Extinction to the nuclear regions of NGC 253 using a foreground dust screen model. 



Ratio: flux(l)/flux(2) 
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Note. — Apertures are rectangular, and centered at the peak of the 2.2 (im continuum brightness (center) and at the Paa source located 
NE of the continuum peak. The errors in the derived extinctions take into account 10% and 20% uncertainties for the nucleus and NE source 
measurements, respectively. 
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the galaxy by Forbes et al. (1993). Other authors (e.g., 
Kcto et al. 1999) claim that a bright non-thermal radio 
source located NE of the peak of the infrared continuum 
emission (source 2 in Turner & Ho 1985) may be the true 
nucleus of NGC 253. It is shown in these figures as a filled 
star symbol. 

The NICMOS Paa line map (Figure 6b) shows very 
bright emission stemming from the peak of infrared con- 
tinuum emission, together with more diffuse emission in a 
plume-like morphology along the plane of the galaxy, pri- 
marily to the NE of the nucleus. As for M82, there is an 
extended component of the Paa emission perpendicular to 
the disk, better seen in our NIC3 images (not shown here) 
which are more sensitive to diffuse emission. This com- 
ponent is however not as strong as that in M82. There 
is a remarkable resemblance of the central Paa emission 
with both the [Neil] 12.8 /im emission line and the nearby 
continuum images presented in Boker, Krabbe, & Storey 
(1998) and Keto et al. (1999). 

4. EXTINCTION 

Both M82 and NGC 253 are known to be highly ob- 
scured (see for instance the comparison between the opti- 
cal HST/WFPC2 F555W and the fiST/NICMOS F222M 
images of M82 in Alonso-Herrero et al. 2001). To calibrate 
the supernova rate in terms of the [Fe n] 1.644 /im emission 
and to obtain the total number of ionizing photons, it is 
necessary to correct the observed quantities for extinction. 

4.1. M82 

We can estimate the extinction to the gas in M82 us- 
ing observations of the ratios of hydrogen recombination 
and near-infrared [Fe n] lines. We used published values 
of the line fluxes for the smallest apertures and matched 
those apertures on the NICMOS line images when mea- 
suring the Paa and [Fen] 1.644 /im fluxes. The intrinsic 
ratios of hydrogen recombination lines are from Hummer 
& Storey (1987) for an electron temperature of T c = 10 4 K 
and electron density of n c = 10 4 cm~ 3 . The ratio between 
[Fell] 1.257 |Um and [Fen] 1.644 /im does not depend upon 
the physical conditions of the gas. We used a polynomial 
fit to the near-infrared extinction law given in He et al. 
(1995). 

Since we do not have NICMOS images of two hydro- 
gen recombination lines, to measure extinction at the full 
NICMOS angular resolution we made use of the broad- 
band filters NIC2 F160W and NIC2 F222M to construct 
a map of extinction to the stars in M82. We assumed a 
color for the unreddened stellar population of H — K = 0.2 
(which applies almost independently of the age of the pop- 
ulation), and, for computational convenience, a foreground 
dust screen model following the extinction law of He et al. 
(1995). The contribution of dust emission at the if -band 
was assumed to be negligible. Figure 7a is the extinction 
map (Ak) from the H — K color map of M82, with the 
positions of the radio sources superimposed. The values 
of the extinction in the if-band derived from this map are 
between Ak = 0.3 mag and Ak = 2.5 mag. The M82 ex- 
tinction map qualitatively looks very similar to those of 
Satyapal et al. (1995), but quantitatively the values of the 
extinction are different. 

The estimates of the extinction to the gas for M82 
are presented in the first part of Table 2 along with the 



line fluxes. The errors in Table 2 arc the uncertain- 
ties associated with the aperture size, aperture centering 
and flux calibration. For comparison, we also estimated 
the extinction through a small aperture (3") using the 
F222M/F160W continuum flux ratio, which gives an es- 
timate of the visual extinction to the stars. As can be 
seen from this table the agreement between the extinc- 
tions to the gas and to the stars is excellent at infrared 
wavelengths. Satyapal et al. (1995) reached a similar con- 
clusion. A much lower level of extinction is indicated by 
comparing the Ha and Paa line fluxes. This behavior is 
known to arise from the different optical depths probed by 
optical and near-infrared lines. This behavior is confirmed 
by the much higher extinction of Ay ~ 27 mag inferred by 
Puxlcy et al. (1989) from the radio H53a to mid-infrared 
Bra line ratio. 

We also computed the extinctions to the El, Wl and 
W2 regions in Achtermann & Lacy (1995) notation (see 
Table 2 for the positions) . These were found to be bright in 
[Neil] 12.8 /im and Bra. We used their Bra measurements 
along with the mid-infrared extinction curve in Rieke & 
Lebofsky (1985). Since Achtermann & Lacy (1995) did 
not quote aperture sizes, we assumed that their fluxes were 
measured over the extent of these regions (approximately 
5" — 6" in diameter). The values of the extinction for 
these three regions are slightly larger than the values we 
would measure from the H — K extinction map, perhaps 
indicating minor differences in the optical depth. 

A further test of the extinction to the gas was based 
on the total number of ionizing photons from the extinc- 
tion corrected Paa line emission map. The number of 
ionizing photons obtained from the extinction-free H53a 
radio recombination line is N-^y = 1 x 10 54 s _1 (Pux- 
ley et al. 1989), and from the 3.3 mm continuum is 
N hy = l.l x lO^s- 1 (Carlstrom & Kronberg 1991). We 
measured through a large aperture the Paa flux on the ob- 
served and the extinction-corrected images. We obtained 
iV Ly = 4 x lO'^s- 1 (observed) and A Ly = 1.4 x lO^s" 1 
(corrected for extinction), for electron temperature T c = 
10 4 K. The total number of ionizing photons corrected for 
extinction is in excellent agreement with a number of esti- 
mates including those of Puxley et al. (1989), Carlstrom 
& Kronberg (1991) and Forster Schreiber et al. (2001). 

4.2. NGC 253 

As for M82, the extinction to the stars was computed 
from the H - K color map of NGC 253 (Figure 7b). It 
shows values between Ay = 3 mag and Ay = 21 mag and 
is very similar to that of Sams et al. (1994) derived from a 
J — K color map. We find a similar value of the extinction 
{Ay ~ 21 mag) for the point with the highest value of the 
extinction which was identified by Sams et al. (1994) as 
the true nucleus of NGC 253. 

To calculate the extinction to the gas in NGC 253 we 
used Paa and Br7 fluxes (the latter from Puxley & Brand 
1995 and Harrison et al. 1998). Both works centered their 
apertures on the approximate coordinates of the secondary 
near-infrared peak and therefore we will assume that the 
largest Br7 fluxes in these works correspond to the primary 
near-infrared peak. If their (0,0) position were to be iden- 
tified with the peak of the continuum emission, their fluxes 
would not be consistent with Forbes et al. (1993) Br7 line 
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emission map. This map and our Paa line emission map 
(Figure 6b) show that the brightest hydrogen line emission 
arises from the peak of the continuum emission. In Table 3 
we give estimates for the extinction for the primary near- 
infrared peak (center) and secondary near-infrared peak 
(positions 3" and 3.9" to the NE of the primary peak for 
Puxley & Brand 1995 and Harrison et al. 1998 works re- 
spectively). The errors take into account 10% and 20% 
combined flux calibration/centering uncertainties for the 
nucleus and NE source, respectively. 

The values of the gas extinction to these two regions 
as derived from the Paa/Br7 line ratio are consistent (to 
within the errors) with the extinction to the stars from 
the H — K color map. Despite this apparent good agree- 
ment we caution the reader of the aforementioned cen- 
tering uncertainties. We also measured the observed and 
corrected for extinction Paa fluxes which yielded a number 
of ionizing photons of A^y = 2.5 x 10 52 s _1 (observed) and 
iVLy = 1.4xl0 53 s _1 (corrected for extinction), for the field 
of view shown in Figure 6b. The latter value is consistent 
with Engelbracht et al. (1998) and Keto et al. (1999) es- 
timates of the number of ionzing photons from their large 
aperture Br7 and [Ne n] measurements, but smaller than 
Puxley et al. (1997) estimate (N Ly = 3.7 x 10 53 s" 1 ) from 
millimeter-wavelength hydrogen recombination lines. 

5. THE [FEII]1.644/nM EMISSION OF M82 AND NGC 253 

5.1. The Distribution of the [Fe Il/l. 644 /jm and radio 

SNRs 

To determine whether there is a spatial correspondence 
between the compact radio sources and the [Fe n] sources 
in both M82 and NGC 253, very precise coordinates are 
needed. Unfortunately, due to problems with the HST 
telemetry during the observations of M82, and the un- 
certainties when producing and rotating the mosaics, the 
coordinates of the NIC3 images may not be completely 
reliable. In addition, it is expected that there are minor 
discrepancies between the radio and the optical/infrared 
coordinate reference systems. 

Our approach was to determine the coordinates of the 
[Fe n] image by trying to find the highest possible number 
of coincidences between the positions of the radio SNRs 
and the [Fen] sources. We use DAOFIND within the 
DAOPHOT package in iraf to find compact sources. The 
term "compact" is not used in the sense of point sources. 
In fact, most of these sources seem to be slightly resolved, 
although due to the under-sampling of the NIC3 PSF, 
some uncertainties remain. The sources were selected to 
be easily distinguishable from the extended diffuse emis- 
sion detected above and below the plane of the galaxy. 

In M82 we identified 65 "compact" sources in the 
[Fe n] 1.644 /im map (not corrected for extinction). For the 
radio sources in M82 we used the positions of SNRs iden- 
tified by searching for compact features within radio maps 
at 408 MHz, 1.4 GHz, 5 GHz and 8.4 GHz, with peak flux 
densities above 1 mJy (Dr. Karen Wills, private commu- 
nication 2000). We excluded from this list radio sources 
that are classified as H n regions based on their radio spec- 
tral index (see McDonald ct al. 2002 for details), resulting 
in 44 radio SNRs. We began the comparison of positions 

1 Since the photometry was done on the [Fell] image registered to 
detected sources (71) is slightly different from that in Section 5.1 



using the coordinates in the header of the NICMOS NIC3 
F164N images and made small offsets of the [Fe n] image, 
cross-correlating the positions of the [Fell] sources with 
the radio SNRs. Based on the resolution of the NIC3 
images and the average angular size of the radio SNRs 
(~ 0.2", Huang et al. 1994), we counted a coincidence 
when the angular separation between a radio source and 
an [Fell] source was < 0.5". With a small shift, we found 
21 coincidences between radio SNRs and [Fen] sources, 
whereas if we used the original coordinates we could only 
find 7 coincidences. Monte Carlo simulations showed that 
the probability of finding 21 coincidences by pure chance 
is 2 x 10~ 6 (see Whitmore & Zhang 2002 for details on 
the method), so we believe that now the infrared and the 
radio images are properly aligned. 

For NGC 253 we found discrepancies between the coor- 
dinates derived for the near-infrared peak of emission from 
the NIC2 and NIC3 images (which were taken during two 
different orbits). Although no problems were reported for 
these two sets of observations, some error must be affect- 
ing the coordinates given in the headers of the images. 
Again as with M82, we cross-correlated the positions of 
the [Fen] 1.644 /xm (36 sources) and 35 radio sources taken 
from Ulvestad & Antonucci 1991 — excluding those sources 
that are classified as H II regions in Ulvestad & Antonucci 
(1997) — using the coordinates of the NIC3 images. We 
find a total of 9 coincidences for a separation between the 
radio sources and [Fell] sources < 0.5". The probability 
of finding 9 coincidences in by pure chance in NGC 253 is 
approximately 9 x 10~ 3 . 

One important conclusion, despite the uncertainties of 
the image alignment, is that we cannot assign an [Fen] 
source to each radio SNR or vice versa. We find approxi- 
mately 30 - 50% of radio SNRs in M82 and NGC 253 have 
an [Fe ii] . This point is illustrated in the contour images of 
the [Fen] 1.644 /xm line emission and continuum emission 
at 1.66 /an in Figures 2a and 2b for M82 and Figures 5a 
and 5b for NGC 253. In all figures we show the positions 
of the radio SNRs as asterisks. 

We also tried to determine if there is a relation between 
the presence of a coincidence and the age of the radio SNR. 
Such a correlation would be expected if the [Fe ii] emission 
were brighter in older supernova remnants. For M82 we 
used the angular sizes of the radio SNR as an indication 
of their age (see Huang et al. 1994). In general, we did 
not find that older radio SNRs tend to have a coincidence 
with an [Fen] source. 

5.2. Luminosity and Origin of the /Fell/1.644 \xm sources 

5.2.1. Supernova Excitation 

We performed aperture photometry on the [Fe ii] sources 
identified in M82 1 and NGC 253 using a fixed 0.8" diame- 
ter aperture. This aperture corresponds to a physical size 
of 12 pc and 10 pc for M82 and NGC 253, respectively. 
Figures 9a and 9b show the distributions of the measured 
[Fen] 1.644 /jm luminosities (not corrected for extinction). 
Some of the sources may be larger than the aperture used 
for the photometry; on the other hand, patchy extinction 
will make some of the sources appear smaller. From in- 
spection of typical source sizes, however, the luminosities 

match the field of view and position of the Pact image, the number of 
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would be accurate on average. Recent studies of Galactic 
SNR show that the [Fe n] emission is filamentary and have 
sizes of the order of 11 pc (e.g, IC443 Rho et al. 2001; see 
also Morel et al. 2002 for SNR in M33). 

For the assumed distances to M82 and NGC 253 the ob- 
served [Fen] 1.644 /im luminosities of these sources range 
from 2.2 x 10 27 W to 1.2 x 10 30 W. As a consistency check 
for the photometry, we measured the [Fen] 1.644 /im fluxes 
for the Fel, Fe2 and Fe3 sources in M82 identified by 
Greenhouse et al. (1997) using the aperture diameters 
listed in their table 1. The agreement with Greenhouse et 
al. (1997) fluxes for these sources is excellent, to within 
3%. However, when we compared our photometry for the 
sources reported in Forbes et al. (1993) for NGC 253 we 
found some discrepancies. 

The corrected [Fen] 1.644 /xm luminosities for compact 
sources in M82 and NGC 253 range from 3 x 10 29 W 
to 2.2 x 10 31 W. As a comparison, Galactic SNRs yield 
[Fen] 1.644 /xm luminosities of 1 x 10 26 W to 7 x 10 28 W 
(Oliva et al. 1989; Keller et al. 1995), SNRs in the LMC 
yield luminosities from 1 x 10 28 W to 3 x 10 29 W (Oliva 
et al. 1989) and SNRs in M33 (assuming a distance of 
720 kpc) 5 x 10 27 to 3 x 10 29 W (Lumsden & Puxley 1995 
and Morel et al. 2002). 

The individual sources account for some 22% and 27% 
of the total [Fen] 1.644 /im emission (extinction corrected) 
measured in M82 and NGC 253, respectively. Given the 
estimates that the [Fen] emitting phase of a SNR has a 
duration of ~ 10 4 yr (Lumsden & Puxley 1995; Morel et al. 
2002) and the supernova rates in M82 and NGC 253 are 
~ 0.10 yr -1 (Huang et al. 1994 and Ulvestad & Antonucci 
1997), the detected [Fell] sources account for ~ 1000 yr 
worth of supernovae. Presumably much of the remaining 
unresolved [Fen] arises from additional supernovae that 
are packed so densely that they are confused in our im- 
ages. Also, SNRs in these two galaxies may expand and 
merge into a general hot ISM within a few thousand years 
losing their identity as individual sources. Photometry on 
the [Fe n] emission arising from the disk of M82 shows that 
as much as 70% of the total emission (disk plus component 
perpendicular to the disk) is associated with SNRs. This 
together with the good agreement between the radio and 
[Fe n] morphologies of the components perpendicular to 
the disk gives further indication for a SNR origin of the 
[Fe II] emission in starburst galaxies. 

Lumsden & Puxley (1995) pointed out that one possi- 
bility is that the [Fe n] lines reach maximum luminosity 
when the remnant as a whole undergoes the transition 
from adiabatic expansion to radiative expansion. Accord- 
ing to Shull & Draine (1986), this happens when radiative 
cooling occurs faster than dynamical times, around 10 4 yr 
after the supernova explosion. Therefore if the brightest 
[Fell] emission occurs at approximately 10 4 yr, it is not 
surprising that we do not find a correspondence for all 
the [Fen] sources and all the radio SNRs. Muxlow et al. 
(1994) found that half of the radio sources in M82 arc un- 
resolved indicative of a young age, but the other resolved 
half are probably intermediate age. This together with the 
fact that we found that 30 — 50% of the radio SNR have 
an [Fen] counterpart implies that the [Fen] sources may 
be^l0 4 yr. 

5.2.2. Excitation in Hu Regions 



One possibility for the lack of one-to-one coincidences 
between the radio SNRs and the [Fe n] sources would be 
that most of the [Fen] emission in these two galaxies is 
produced in Hll regions (photoionization). To check this, 
we measured the Paa fluxes for the sources detected in 
the [Fen] extinction corrected map of M82, and assumed 
an [Fen] 1.644 /im/Paa ratio for Hu regions as that ob- 
served in Orion. We found that only ~ 6 — 8% of the 
total [Fen] 1.644 /im arises from Hu regions. This further 
supports the hypothesis that most of the [Fe n] 1.644 /xm 
emission along the plane of M82 is produced by SNRs. 

The [Fe Ii]/Br7 ratios in H n regions and SNR differ by a 
factor of a hundred or more, and therefore can be used to 
indicate the excitation mechanism of the [Fe n] line within 
our images. In Alonso-Herrero et al. (1997) we showed 
that the behavior of the [Fe n]1.644/im/Br7 line ratio in 
starbursts and active galaxies can be understood as a pro- 
gression from pure photoionization to pure shock excita- 
tion going from pure H n regions through starbursts, tran- 
sition objects, Seyferts to supernova remnants. 

Figures 8a and 8b are the observed [Fen] 1.644 /um/Paa 
maps for M82 and NGC 253, respectively. The values 
of the [Fe n]1.644/im/Paa ratio range from 0.02 in black 
(mostly along the disk of the galaxy) to 0.8 in white (some 
of the bright [Fen] sources). These ratios translate into 
[Fe n] 1.644 /im/Br7 line ratios of between 0.24 and 10 (see 
also the histograms in Figure 10). Since these maps are 
not corrected for extinction, these numbers are actually 
lower limits. When the positions of radio sources are plot- 
ted in the [Fen] 1.644 /im/Paa maps we find again that 
there is a number of coincidences between the radio SNRs 
and the regions of enhanced [Fen] 1.644 /im emission. This 
indicates that the relative contribution of the SNRs is 
nearly as large in the diffuse background as it is in the 
radio SNR. Interestingly, we find that most of the coinci- 
dences with the radio SNRs occur in regions of enhanced 
[Fe n] 1.644 /im/Paa line ratio, or conversely we do not find 
a large number of coincidences in regions with high Paa 
flux (that is, Hu regions). 

6. CALIBRATION OF THE SUPERNOVA RATE - [FEII] 
LUMINOSITY RELATION 

The supernova rate in M82 has been estimated from 
number counts of the radio supernova remnants and their 
ages to be 0.11 ± 0.05 yr" 1 (Huang et al. 1994). The 
total [Fell] 1.644 /im luminosity in M82 is L([Feil]) to tai = 
1.9 x 10 33 W. Our value for the [Fen] luminosity of M82 
is approximately 1.5 and 2 times higher than the previous 
estimates by Greenhouse et al. (1997) and Vanzi & Rickc 
(1997), respectively. This difference is most likely due to 
better extinction estimates made possible by our NICMOS 
images. If the supernova rate in M82 is 0.11 yr -1 , then we 
can derive the following calibration for the supernova rate 
in terms of the [Fe II] 1.644 /xm luminosity: 



(SNr) 



M82 



0.6 x 



L 



[Foil] 



(yr" 1 ) 



(1) 



10 34 W 

The total [Fen] 1.644 /im luminosity (corrected for extinc- 
tion) measured in NGC 253 is 6.1 x 10 32 W. An upper 
limit on the supernova rate from radio supernova counts 
is < 0.1 - 0.3 yr -1 (Ulvestad & Antonucci 1997). A rate 
of 0.08 yr" 1 has been calculated assuming that the non- 
thermal radio flux is entirely due to synchrotron emission 



8 



The [Fe n] emission in M82 and NGC 253 



(Van Buren & Greenhouse 1994); this value is consistent 
with the similar number of compact radio sources in this 
galaxy compared with M82. The calibration of supernova 
rate in terms of [Fell] 1.644 /xm luminosity for this galaxy 
is then 



(SNr) 



NGC253 



1.3 x 



L 



[Foil] 



(yr- 1 ) 



(2) 



10 34 W 

Since M82 is the more thoroughly studied galaxy in the 
radio, we give higher weight to the calibration for it in 
averaging these results to estimate 



(SNr) = 0.8 x 



J [Fon] 



(y^ 1 ) 



(3) 



10 34 W 

In the above calculations, we have included the entire 
[Fe n] luminosity, since that is the parameter usually mea- 
sured in a distant galaxy. However, in these nearby ex- 
amples, we can compare the [Fen] emission produced in 
H II regions with that from supernovae. We assume that 
the ratio for a typical Hn region is /([Fen] 1.644 /xm)Hii — 
6.6 x 10~ 3 .f(Paa) (e.g., Orion). We obtained that the 
[Fe n] luminosity produced in H n regions accounts for less 
than 8% the total [Fell] luminosity in M82. 

7. SUMMARY AND CONCLUSIONS 

We have presented ffST/NICMOS observations of the 
starburst galaxies M82 and NGC 253. The morphology of 
the [Fell] 1.644 /xm line emission in both galaxies displays 
two components: (1) a disk component with a number of 
superimposed compact sources and (2) a diffuse compo- 
nent in the form of filaments above and below the disk of 
the galaxy. In M82 we find a remarkable similarity be- 
tween the extended component, the fingers or chimneys 
above and below the disk, of the radio and [Fen] emis- 
sions. 

In M82, we have detected in both Paa and [Fell] a nu- 
clear ring of star formation. This ring, with an approxi- 
mate diameter of 400 pc, has also been observed in other 
near- and mid-infrared hydrogen recombination lines and 
in the [Neii]12.8/xm line. The peak of the 2.2 /xm contin- 
uum brightness of M82 lies very close to (within 1") the 
dynamical center, and corresponds to the center of the ring 
of star formation. In NGC 253, bright Paa emission stems 
from the near-infrared continuum peak and extends in a 
plume-like morphology along the disk of the galaxy. There 
is also a diffuse Paa component above and below the plane 
of galaxy, however, much fainter than that in M82. 

Prior to calibrating the surpernova rate in terms of the 
[Fe n] emission, we have derived the extinction to the cen- 



tral regions of NGC 253 and M82. To assess the extinction 
correction, we have compared the number of ionizing pho- 
tons infered from the extinction-corrected Paa maps with 
other estimates based on longer wavelength indicators, and 
found a good agreement. 

Since the [Fe n] emission in starburst galaxies is thought 
to be associated with supernova activity, we have cross- 
correlated the positions (with a spatial resolution of a few 
tenths of arcseconds) of the radio SNRs and the identified 
compact [Fen] sources in M82 and NGC 253. We have 
found that approximately 30 — 50% of the radio SNRs 
have an [Fen] counterpart. However, we cannot associate 
an [Fell] source to each radio SNR or vice versa. This is 
well understood in terms of the different ages traced by 
the radio SNRs (a few hundred years) and the [Fe n] stage 
of SNRs (<J 10 4 yr). The compact [Fen] emitting sources 
in M82 and NGC 253 only contribute to some 20 - 30% 
of the total [Fen] 1.644 /im emission. However, much of 
the remaining emission in the disk of the two galaxies is 
most likely produced by SNRs that expanded and merged 
into a general ISM a few 10 4 yr ago. Photometry on the 
[Fe n] 1.644 /im emission arising from the disk of M82 shows 
that as much as 70% of the total [Fe n] 1.644 /im emission 
(corrected for extinction) is associated with SNRs. The 
remaining diffuse [Fe n] emission perpendicular to the disk 
of M82 seems to be related to the superwind (via expulsion 
of material from the central starburst), as argued from the 
radio emission morphology (Wills et al. 1999). 

All these arguments support the supernova origin for 
the observed near- infrared [Fe n] line emission in starburst 
galaxies. Thus we have recalibrated the supernova rate in 
terms of the [Fe n] luminosity relation for starburst galax- 
ies using our newly measured (corrected for extinction) 
[Fen] 1.644 /xm luminosities of M82 and NGC 253. 
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Fig. 1a. — M82: continuum emission at 1.66 /im (NIC3 F166N image), displayed on a logarithmic scale. 
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Fig. 1b. — M82: [Fell] 1.644 /im line emission (continuum subtracted NIC3 F164N image), displayed on a logarithmic scale. 



Fig. 2a. — Contour map of the M82 continuum emission at 1.66 fim on a linear scale. The asterisks are the positions of the radio SNRs 
detected in radio maps (Karen Wills, private communication 2000, see text for details). 
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Fig. 2b. — Contour map of the M82 [Fe II] 1.644 fim line emission (continuum subtracted) on a linear scale. Symbols as in Figure 

figure3b.jpeg 

Fig. 3a. — M82: continuum emission at 1.90 (im (NIC2 F190N images) displayed on a logarithmic scale. 

figure3b.jpeg 

Fig. 3b. — M82: Paa (1.87 ^m) line emission map (continuum subtracted NIC2 F187N images) displayed on a logarithmic scale. 

figure4a.jpeg 

Fig. 4a. — NGC 253: continuum emission at 1.66 fim (NIC3 F166N image), displayed on a logarithmic scale. 
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figure4b.jpeg 



Fig. 4b. — NGC 253: [Fell] 1.644 /im line emission (continuum subtracted NIC3 F164N image), displayed on a logarithmic scale. 




Fig. 5a. — Contour map of the NGC 253 continuum emission at 1.66 ^tm on a linear scale. The asterisks are the positions of the 6 cm 
radio SNRs from Ulvestad & Antonucci (1991) with S pea k > 1 mjy, excluding those sources classified as Hll regions in Ulvestad & Antonucci 
(1997). The filled triangle is the position of the peak of the near-infrared emission, whereas the filled star symbol is the position of the 
non-thermal radio source (see text for details). 
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Fig. 5b. — Contour map of the NGC 253 [Fe II] 1.644 (im line emission (continuum subtracted) on a linear scale. Symbols as in Figure 5a. 
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Fig. 6a. — NGC 253: continuum emission at 1.90 fim (NIC2 F190N images) displayed on a logarithmic scale. 
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Fig. 6b. — NGC 253: Paa (1.87 //m) line emission map (continuum subtracted NIC2 F187N images) displayed on a logarithmic scale. 
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Fig. 7a. — M82: H - K color map (NIC2 F160W - NIC2 F222M) which represents the extinction to the stars. Symbols arc as in 
Figure 2a. Dark regions indicate regions of higher extinction. The circular black spot on the upper right corners of the individual images is 
the coronographic hole. 
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Fig. 7b.— NGC 253: H — K color map (NIC2 F160W - NIC2 F222M). Symbols are as in Figure 5a. The black spot on left bottom corners 
of the individual images is the coronographic hole. 
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Fig. 8a. — M82: [Fen]1.644/im/Paa line ratio map. Lighter colors indicate regions of enhanced [Fell]1.644/im/Paa line ratios. Superim- 
posed are the positions of the radio SNRs as in Figure 2a. 
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Fig. 8b. — NGC 253 [Fell] 1.644 /^m/Paa line ratio map. Colors as in Figure 8a. Superimposed are the positions of the radio SNRs as in 
Figure 5a. 
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Fig. 9. — Distributions of the [Fell] 1.644 (im luminosities (not corrected for extinction) of the compact sources in M82 and NGC 253 
measured through a 0.8" diameter aperture. 
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Fig. 10. — Histograms of the observed [Fell] 1.644 fim/Hr-y line ratios of the compact [Fell] sources. 
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